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Abstract 

In this paper we demonstrate that neutral Dirac particles in external electric fields, 
which are equivalent to generalized Dirac oscillators, are physical examples of quasi- 
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of the system based on the sl(2) symmetry are discussed separately in spherical, cylin- 
drical, and Cartesian coordinates. Some exactly solvable field configurations are also 
exhibited. 
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I. Introduction 

Search for exact solutions of wave equations, whether non relativistic or relativistic, has 
been there since the birth of quantum mechanics. In the case of non relativistic wave equa- 
tions, i.e. Schrodinger equations, the class of potentials for which the complete spectrum or a 
part of the spectrum are known has grown over the years. However in the case of relativistic 
wave equations, in particular, the Dirac equation very few exactly solvable electromagnetic 
field configurations like homogeneous magnetic fields [1], homogeneous electrostatic fields 
[2], constant parallel magnetic fields [3] etc. are known. 

It may be mentioned that while exact solvability is desirable, in practice it is not always 
possible to determine the whole spectrum. In non relativistic quantum mechanics a new 
class of potentials which are intermediate to exactly solvable ones and non solvable ones have 
been found recently. These are called quasi-exactly solvable (QES) problems for which it is 
possible to determine algebraically a part of the spectrum but not the whole spectrum [[4]- 
[9]]. Although there are exceptions, usually a QES problem admits a certain underlying Lie 
algebraic symmetry which is responsible for the quasi-exact solutions. For Dirac equations, 
it had been shown that Dirac equation with Coulomb interaction supplemented by a linear 
radial potential [10], and planar Dirac equation with Coulomb and homogeneous magnetic 
fields [11] are QES systems. More recently, the quasi-exact solvability of the Pauli equation 
was exhibited [12]. This immediately implies that the Dirac equation coupled minimally to a 
vector potential is QES, since the square of the Dirac Hamiltonian in this form is proportional 
to the Pauli equation, up to an additive constant, namely, the square of the rest energy. 

Interactions of spin- 1/2 fermions with electromagnetic fields are usually introduced in 
the Dirac equation following the minimal coupling prescription. However, this prescription 
is not the only possibility. For instance, the interaction due to anomalous magnetic moment 
is incorporated via the non-minimal coupling procedure, the resulting equation being called 
the Dirac-Pauli equation [13]. An interesting example is neutral fermions interacting with 
electromagnetic fields. There are a number of interesting phenomena involving non mini- 
mal coupling of neutral fermions to external electric fields. We recall that, the celebrated 
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Aharonov-Casher effect [14], which has been observed experimentally [15] is an example 
of such an interaction. Another interesting system involving neutral fermions is the Dirac 
oscillator, which is of considerable interest in quantum chromodynamics [16]. However, un- 
like the case with minimal coupling, studies in exact solutions of the Dirac equations with 
non-minimal couplings are rather scanty [17, 18, 19], not to mention studies in quasi-exact 
solutions of these systems. Only recently, efforts have been directed to exploring certain 
structure relating to such systems. For instance, the underlying supersymmetry (SUSY) of 
the Aharonov-Casher system [20], as well as that of a system of neutral fermions interacting 
with a central electric field were pointed out [21]. 

In the present paper we shall examine quasi-exact solvability of Dirac equation coupled 
non-minimally to external electric fields. In particular, based on the underlying s/(2) sym- 
metry we shall determine the forms of electric fields which give rise to QES Dirac-Pauli 
equations. 

The organization of the paper is as follows. In Sect. II the Dirac-Pauli equation is pre- 
sented, and its equivalence with the generalized Dirac oscillators demonstrated. Sect. Ill dis- 
cusses the exact and quasi-exact solvability of the Dirac-Pauli equation with electric fields 
in the spherical coordinates. Sect. IV and V are devoted to the same problems with electric 
fields in the cylindrical and Cartesian coordinates, respectively. Sect. VI concludes the paper. 

II. The Dirac-Pauli Equation for neutral Dirac particles 

We consider the motion of a neutral fermion of spin-1/2 with mass m and an anomalous 
magnetic moment /z, in an external electromagnetic field described by the field strength 
F^y. The fermion is described by a four-component spinor \1/ which obeys the Dirac-Pauli 
equation [13] 



(17"^ - -fJur^F^ - m)tt = , 



(1) 



where 7 M = (7°,7) are the Dirac matrices satisfying 



(2) 
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with g» v = diag(l, -1, -1, -1), and 



(3) 



In terms of the external electric field E and magnetic field B, one can rewrite the second 
term in eq.(l) as 



—a^Ffu, = ia E - £ B , 



(4) 



where a. = 7°7, and S fe = ^e kl ^a lJ . Here e fcy is the totally anti-symmetric tensors with 
e 123 = I. For time-independent fields one may set 



tf(t,r) = e- i£ V(r) , 

and eq.(l) becomes 

Hi( = Sip, 
with the Hamiltonian H being given by 

H = a ■ p + i/ij ■ E - /!/?£ ■ B + (3m , 



(5) 



(6) 



(7) 



where p = — iV and /3 = 7 . 

In this paper we shall consider cases with only electric fields described by the the Hamil- 
tonian 



H = ct ■ p + ifij ■ E + (3m . 
We choose the Dirac matrices in the standard representation 



(8) 
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(9) 

where cr are the Pauli matrices. We also define ip = (%, </?)*, where t denotes transpose, and 
both x an d V 9 are two-component spinors. Then the Dirac-Pauli equation becomes 



cr • (p - i/iE)x = (8 + m)v? , 
cr • (p + i/iE)v? = (8 - m)x, 
4 



(10) 



In the next three sections, we will demonstrate the exact and quasi-exact solvability of 
this set of equations in the spherical, cylindrical, and Cartesian coordinates. 



Relation with the Dirac oscillator 



Here we show that the above system is equivalent to the generalized Dirac oscillators. 
Since its introduction by Moshinsky and Szczepaniak in 1989, the original Dirac oscillator has 
attracted much attention in recent years [16]. It was introduced by replacing the momentum 
p in the field-free Dirac equation by p — imuj3r, where m is the mass of the particle and u 
is the oscillation frequency. The Hamiltonian of this system is thus given by 



This system is exactly solvable. 

Now if one compares eq.(ll) with (8), one sees immediately that the Dirac oscillator is 
simply a special case of neutral Dirac particle in external electric field, if one makes the 
transformation 



That is, when the electric field E is proportional to the position vector r. Hence, with a 
general r-dependent electric field E(r), the system described by the Hamiltonian (8) can be 
viewed as generalizations of the original Dirac oscillator. In the discussions below, we shall 
concentrate on the neutral Dirac particles. But all conclusions thus obtained are applicable 
directly to the generalized Dirac oscillators. 

III. Electric fields in spherical coordinates 

First let us consider central electric field E = E r r. In this case, one can choose a complete 
set of observables to be (H, J 2 , J 2 ,S 2 = 3/4, K). Here J is the total angular momentum 
J = L + S, where L is the orbital angular momentum, and S = |S is the spin operator. 
The operator K is defined as K = 7°(S • L + 1), which commutes with both H and J. 



H = ot • (p — imu/3r) + f3m . 



(11) 



yuE — > mur . 



(12) 
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Explicitly, we have 



K = diag (k, —kj , 
k = a ■ L+ 1 . 



The common eigenstates can be written as [22] 



r 



i f / 



here yf m .(0,(f)) are the spin harmonics satisfying 



J 2 y k 



1 3 



jmj J \J ' "J^jrrij ' J 2' 2' 



= 77? -"V* 



fc = ±0" + \) > 



and 



where r is the unit radial vector. Using the identity 

er • p = i(cr • r) ^-d r + ^(cr ■ L)^ , 

one gets 

cr • (p ± i/iE) = i(cr • r) (^-d r + — ± liE r ^j . 

With these, eq.(10) reduces to 



U + - + ^)f- 
\ar r J 

—^ + - + piE r ) A 
ar r J 



(£ + m)f + , 
(£-m)/_ . 



(13) 



(14) 



(15) 
(16) 
(17) 



(18) 



(19) 



(20) 



(21) 
(22) 



This shows that /_ and / + forms a one-dimensional SUSY pairs [21]. For if we write the 
superpotential W as 



W = - + fiE r , 
r 



(23) 



then eqs.(21) and (22) become 



with 



A'A + f_ = (S 2 -m 2 )f_, (24) 
A + A~U = (S 2 -m")U, (25) 



A ± = ±4- + W . (26) 
dr 

Explicitly, the above equations read 

(-^ + W 2 TW^U=(S*-m>)f T . (27) 

Here and below the prime means differentiation with respect to the basic variable. Eq.(27) 
clearly exhibits the SUSY structure of the system. The operators acting on f± in eq.(27) 
are said to be factorizable, i.e. as products of A~ and A + . The ground state, with £ 2 = m 2 , 
is given by one of the following two sets of equations: 

A + ^\r) = , /| 0) (r)=0; (28) 
A-ff\r) = , / W (r)=0, (29) 

depending on which solution is normalizable. The solutions are generally given by 

f T oc r Tk exp (t J drfiE^j . (30) 

We now classify the forms of the electric field E r (r) which allow exact and quasi-exact 
solutions. To be specific, we consider the situation where k < and / drjj,E r > 0, so that 
/I - 1 is normalizable, and = 0. The other situation can be discussed similarly. In this 
case, eq.(23) becomes 

W = -— + fiE r . (31) 
r 

We determine the forms of E r that give exact /quasi-exact energy £ and the corresponding 
function /_. The corresponding function / + is obtained using eq.(21). 
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A. Exactly solvable cases 

Comparing the forms of the superpotential W in eq.(31) with Table (4.1) in [23], one 
concludes that there are three forms of E r giving exact solutions of the problem : 

i) oscillator-like : fJ>E r (r) oc r ; 

ii) Coulomb potential-like : iiE r {r) oc constant ; 
hi) zero field-like : /iE r (r) oc 1/r . 

Case (i) and (ii) had been considered in [18] and [19], and case (hi) in [18]. 
We mention here that the case with oscillator-like field, i.e. case (i), is none other than 
the spherical Dirac oscillator. 

B. Quasi-exactly solvable cases 

The form of the superpotential W in eq.(31) fits into three classes, namely, Classes VII, 
VIII and IX, of s/(2)-based QES systems in [5]. Electric field configurations permitting any 
number of solvable excited states in each of these classes can be constructed according to 
the general procedure given in [12]. The procedure makes use of the connection between 
quasi-exact solvability and the SUSY structure (or, equivalently, the factorizability) of the 
systems. Below we describe the construction of QES E r for the simplest cases in Class VII. 
Other classes can be considered similarly. 

An outline of the method of construction 

Here we describe briefly the main ideas underlying the construction of QES systems as 
given in [12]. As before, we shall concentrate only on solution of the upper component /_, 
which is assumed to have a normalizable zero energy state. 

Eq.(27) shows that /_ satisfies the Schrodinger equation H_f- = e/_, with energy pa- 
rameter e = £ 2 — m 2 , and 

H = A~A + 

= -£z+V(r)> (3 2 ) 
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with 

V(r) = W{rf - W\r) . (33) 

We shall look for V(r) such that the system is QES. According to the theory of QES models, 
one first makes an imaginary gauge transformation on the function /_ 

/_(r) = 4>{r)e-M , (34) 

where g(r) is called the gauge function. The function 0(r) satisfies 

d^(r) + 2gl d_^ + r (r) + gl , _ g/2 , 0(r) = £0(r) _ (35) 



dr 2 dr 

For physical systems which we are interested in, the phase factor exp(— g(r)) is responsible for 
the asymptotic behaviors of the wave function so as to ensure normalizability. The function 
(j){r) satisfies a Schrodinger equation with a gauge transformed Hamiltonian 

H G = -^ + 2W (r)-^+[v(r) + W^-W*} , (36) 

where W (r) = g'(r). Now if V(r) is such that the quantal system is QES, that means 
the gauge transformed Hamiltonian H G can be written as a quadratic combination of the 
generators J a of some Lie algebra with a finite dimensional representation. Within this finite 
dimensional Hilbert space the Hamiltonian Hq can be diagonalized, and therefore a finite 
number of eigenstates are solvable. For one-dimensional QES systems the most general Lie 
algebra is s/(2) ([4]- [9]). Hence if eq.(36) is QES then it can be expressed as 

H G = J2 C ah J a J b + ]T C a J a + constant , (37) 

where C ab , C a are constant coefficients, and the J a are the generators of the Lie algebra 
s/(2) given by 

J+ = z 2 ^--Nz, 

dz 

= *|-£, "=0,1,2... 
J - - I ■ < 38 > 



Here the variables r and z are related by z = h(r), where h(-) is some (explicit or implicit) 
function . The value j = N/2 is called the weight of the differential representation of sl(2) 
algebra, and iV is the degree of the eigenfunctions 0, which are polynomials in a (N + 1)- 
dimensional Hilbert space 

= (z - zi)(z - z 2 ) ■ ■ ■ (z - z N ) . (39) 

The requirement in eq.(37) fixes V(r) and W (r), and H G will have an algebraic sector 
with N + l eigenvalues and eigenfunctions. For definiteness, we shall denote the potential V 
admitting N + 1 QES states by Vn- From eqs.(34) and (39), the function /_ in this sector 
has the general form 

f- = (z- Zl )(z -z 2 )---{z- z N ) exp (- J" W (r)dr^j , (40) 

where z\ (i = 1,2, ...,N) are parameters that can be determined by plugging eq.(40) 
into eq.(35). The algebraic equations so obtained are called the Bethe ansatz equations 
corresponding to the QES problem [9, 11]. Now one can rewrite eq.(40) as 

/- = exp (- f W N (r, {z t })dr^j , (41) 

and 

W N (r, {zi}) — W (r) - ^ TT^ZTT ■ ( 42 ) 

i=i ) z i 

There are N+l possible functions Wn(t, {zi}) for the N+l sets of eigenfunctions 0. Inserting 
eq.(41) into H_f_ = e/_, one sees that Wjy satisfies the Ricatti equation [8, 24] 

W 2 N - W' N = V N - e N , (43) 

where is the energy parameter corresponding to the eigenfunction /_ given in eq.(40) for 
a particular set of iV parameters {zi}. 

From eqs.(32), (33) and (43) it is clear how one should proceed to determine the electric 
fields so that the Dirac-Pauli equation becomes QES based on si (2): one needs only to 
determine the superpotentials W(r) according to eq.(43) from the QES potentials V(r) 
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classified in [5]. This is easily done by observing that the superpotential W corresponding 
to N = is related to the gauge function g(r) associated with a particular class of QES 
potential V(r) by g'(r) = Wo(r). Once Wo is obtained, the corresponding electric field E^ 
is obtained through eq.(31): 

liE^ = W + \k\/r . (44) 

This is the required electric field that allows the weight zero (j — N — 0) state, i.e. the 
ground state, to be known in that class. The more interesting task is to obtain higher weight 
states (i.e. j > 0), which will include excited states. For weight j (N = 2j) states, this 
is achieved by forming the superpotential W N (r, {z,j}) according to eq.(42). Of the N + 1 
possible sets of solutions of the Bethe ansatz equations, the set of roots {z\, z 2 , ■ ■ ■ , z N } to be 
used in eq.(42) is chosen to be the set for which the energy parameter of the corresponding 
state is the lowest (usually it is the ground state). The required electric field which gives 
rise to the N + 1 solvable states is then obtained as 

fiEW = w N + \k\/r . (45) 

For the spherically symmetric electric fields which we consider in this section, there are 
three possible types of s/(2)-based QES field configurations. They belong to Class VII, VIII, 
and IX in the classification listed in [5]. Below we shall illustrate our construction of QES 
electric fields through Class VII QES systems. 

Class VII 

The general potential in Class VII has the form 

V N {r) = a 2 r 6 + 2abr A + [b 2 - a (4N + 2 7 + 3)] r 2 + 7 (7 - 1) r' 2 - b (2 7 + 1) , (46) 
where a, b and 7 are constants. The gauge function is 

g(r) = jr 4 + -r 2 — 7 lnr . (47) 
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We must have 0,7 > to ensure normalizability of the wave function. Eqs.(47) and (31), 
together with the relation W (r) = g'(r), give us the electric field E^: 

/i^ 0) (r) = ar 3 + br . (48) 

The Dirac-Pauli equation with this field configuration admits a QES ground state with 
energy S 2 = m 2 (e = 0) and ground state function /_ oc exp(— go(r)). Also, here we have 
7 = \k\. We retain the symbol 7 so that some general formulae in this section can be carried 
over in the next section simply with 7 redefined. 

To determine electric field configurations admitting QES potentials Vjv with higher 
weight, we need to obtain the Bethe ansatz equations for <f>. Letting z = h(r) = r 2 , eq.(35) 
becomes 



d 2 d 
4z— + (4az 2 + 4bz-2 (2 7 + 1)) — - {4aNz + e) 



<f>(z) = . (49) 



dz 2 ^ ' dz 

In terms of the si (2) generators J + , J~ and J°, the differential operator acting on <f>(z) in 
eq.(49) can be written as 

Ty U = -4J° J' + AaJ + + AbJ° - 2 (N + 2 7 + 1) J~ + constant . (50) 

For N = 0, the value of the e is e = 0. For higher N > and <f>(r) = Y[f=i{z — Zi), the electric 
field E { r N \r) is obtained from eqs.(42) and (45): 

^) W = ^o )(r) _|_^)_. (51) 

For the present case, the roots z^s are found from the Bethe ansatz equations 

2az 2 + 2bz t - (2 7 + 1) - T^77 = ' i = l,-..,N , (52) 
and e in terms of the roots Zj's is 

6 = 2(2 7 + l)Er- ( 53 ) 

For N — 1 the roots Z\ are 



N 1 



i=l Zi 



_ -6±^ + 2a(2 7 + l) 
Zl ~ 2a ' lMj 
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and the values of e are 

e ± = 2 (b ± ^Jb 2 + 2a(2 7 + 1)) . (55) 

For a > 0, the root z{ — — \z[ \ < gives the ground state. With this root, one gets the 
superpotential 

2t t 

W^(r) = ar 3 + 6r-- -— - - . (56) 

r z + \z 1 \ r 

From eq.(51), the corresponding electric field is 

9r 

^«(r) = ar 3 + br- - - . (57) 

r + pi I 

The QES potential appropriate for the problem is 

V(x) = Wl - W[ , 

= Vi-e. (58) 

The one-dimensional SUSY always sets the energy parameter of ground state at e = 0. 
Hence, the ground state and the excited state have energy parameter e = and e = e + — e~ = 
4^/6 2 + 2a(2 7 + l), and wave function 

/_ oc e~ 9o{r) (r 2 - zi) (59) 

and 

/_ oc e~ 9o(r) (r 2 - z+) , (60) 

respectively. 

QES potentials and electric fields for higher degree iV can be constructed in the same 
manner. 



IV. Electric fields in cylindrical coordinates 

We now treat the Dirac-Pauli equation in the cylindrical coordinates. In this case, it is 
easier to solve iJ 2 $ = £ 2 $ instead of Hip = Sip. Once <3> is solved, the desired solution ip is 
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given by ip = (H + E )<E». Setting $ = (/_, /+)*, the equation (H 2 - £ 2 )$ = becomes 

p 2 =F /iV • E + /i 2 £ 2 ±v<t-(2ExV-Vx E)] f T = (S 2 - m 2 ) f T . (61) 
We consider electric field configuration which has cylindrical symmetry, namely, 



E x = xf(p) , E y = yf(p) , E z = , 



(62) 



where p = \Jx 2 + y 2 , and f(p) is some function of p. It can be checked that the conserved 
quantities are p z and J z = L z + S 3 /2. Thus the eigenstates can be chosen as 

/ r> / ^ \ ik z ( 



f- 



R 1 (p)e ll t' 
R 2 (p) e i(W)<t> 



I 



R 3 (p)e u <t> 



(63) 



$ is an eigenstate of p z and J z with eigenvalues k z [k z real) and j z = 1+1/2 (I = 0, ±1, ±2 . . .), 
respectively. Eq.(61) reduce to four decoupled ones: 

p 2 f 2 -2pf(l + l)-ppf> + l -^l\ Rl (p) = cflxG*) 

P 

+ ^p 2 / 2 + 2pfl - + 1 + I 



dp 2 
d 2 



dp 2 
d 2 



+ pp 2 f 2 + 2pf (l + l) +fl pf' + 



l 2 - 

~7 



2 , pp 2 f - 2pfl + ppf + i 



Rs(p) 

R 4 (p) 



eR 4 (p) 



(64) 
(65) 
(66) 
(67) 



dp 2 p* 
where e = £ 2 — m 2 — k 2 . From the above equations it follows that eqs.(64) and (67) represent 
a pair of one- dimensional SUSY partners [12], with superpotential given by 



W(p) = ppf(p) - 1 , 7 = \l\ + 1/2 ■ 
P 



(68) 



Eqs.(65) and (66) form another pair of SUSY partners, obtainable from eqs.(67) and (64), 
respectively, by changing p to —p. Results in [12] can therefore be carried over directly with 
minor modifications. 

To be specific, let us assume pf > 0. For I > 0, we have R 2 = R3 = 0. The ground state 
is given by 



Ri oc exp(-#(p)) , 

/p 
p'f(p')dp-l\np , 



(69) 
(70) 
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and i? 4 = 0. For excited states, the components R\ and i? 4 are related by i? 4 oc A + Ri, where 
A + is defined by eq.(26) and (68). For I < — 1, we simply interchange Ri and R 2 , and i? 3 
and i?4. 

Comparing eq.(68) with (31), it is seen that all results in Sect. Ill can be carried over 
by changing r , E r and 7 = \k\ in eq.(31) to p, pf, and 7 = \l\ + 1/2, respectively. Hence 
we will not discuss this case in details, but only quote some main results. For instance, the 
exactly solvable cases are: 

i) oscillator-like : pf(p) oc constant ; 

ii) Coulomb potential-like : pf(p)ocl/p; 

iii) zero field-like : pf(p)ocl/p 2 . 

As in Sect III, the quasi-exactly solvable cases are Classes VII, VIII and IX of QES 
systems based on the s/(2) algebra as classified by Turbiner [5]. For the case of Class VII, 
the simplest two configurations of the electric fields are obtained from the function f(p) as 
follows: 

i) one-state case: pf^(p) = ap 2 + b (a > 0) ; 

ii) two-state case: pf^\p) = ap 2 + b — 2/(p 2 + \z{ |) , where z{ is given by eq.(54). 

V. Electric fields in Cartesian coordinates 

Finally we consider electric field configurations in one dimension. For definiteness, we 
assume 



The conserved quantities in this case are p y and p z , with eigenvalues k y and k z , respectively. 
Eigenstate can be written as 



E x {x) = W(x) , E y = E z = . 



(71) 





(72) 



\ f+( x ) J 



Substituting eqs.(71) and (72) into (61), we get 



+ (pW(x) T ozk y ) 2 t »W' f T = (S 2 



— m 



2 - u ■ 



(73) 
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Again, these equations are in the same form as eq.(3) in [12] and eq.(408) in [23]. Similar 
to the case in Section IV, the four components of ip are related by two one-dimensional 
SUSY, namely, the 1st (2nd) component of /_ and the 2nd (1st) component of / + form a 
one-dimensional SUSY pair. As such, this system can be dealt with according to the general 
procedures outlined in the last section and in [12] . Here we shall only give the classifications 
of the electric field configurations which permit exact and quasi-exact solutions. 

A. Exactly solvable cases 

The exactly solvable cases have been classified in [23]. There are three types of exactly 
solvable field configurations, namely, 



From the discussions given in [12], the system of eq.(73) are QES based on the si (2) 
algebra. In fact, six classes of QES electric field configurations can be constructed. These 
classes correspond to Class I to Class VI in Turbiner's classification. One can follow the 
general procedures described in Sect. Ill to construct in each class the electric fields which 
admit different numbers of QES states, including excited states. Here we only present the 
main results for Class I, referring the readers to [12] for further details. 



We concentrate on the first set of one-dimensional SUSY system in eq.(73), namely, the 
first component of /_ and the second component of / + . The superpotential for the system 
is 



ax + c (shifted oscillator-like) ; 

c\ exp(— ax) + C2 (Morse-potential); 
tanh(ax) + c (Rosen-Morse II) ; 



(74) 



Here c, ci, C2 and a are constants. 



B. Quasi-exactly solvable cases 



Class I 



W{x) = /iW(x) - k y . 



(75) 
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According to Turbiner's classification, the QES potential belonging to Class I has the 
form 

V N (x) = a 2 e~ 2ax - a [a(2N + 1 ) + 2b] e~ ax + c (2b - a) e ax + c 2 e 2ax + b 2 - 2ac . (76) 

Without loss of generality, we assume a, a, c > for definiteness. The corresponding gauge 
function g(x) is given by 

g ( x ) = - e ~ ax + -e ax + bx . (77) 
a a 

One should always keep in mind that the parameters selected must ensure convergence of 
the function exp(— g(x)) in order to guarantee normalizability of the wave function . The 
potential V(x) that gives the ground state, with energy parameter e = £ 2 — m 2 — k 2 = 0, is 
generated by 

V(x) = V 

= W 2 -W^ (78) 

with 

W (x) = g\x) 

= -ae~ ax + ce ax + b . (79) 

The corresponding electric field is 

fiEi 0) = -aae- ax + cae ax + b + k y . (80) 

To obtain electric fields and the corresponding potentials which admit solvable states 
with higher weights j, we must first derive the Bethe ansatz equations. To this end, let us 
perform the change of variable z = h(x) = exp(— ax). Eq.(35) then becomes 



2 d 2 



2az 2 - {2b + a)z - 2c 



d_ 

dz 



-2aNz - - 
a 



4>{z) = . (81) 



The differential operator acting on <p(z) can be written as a quadratic combination of the 
si (2) generators J + , J~ and J° as 

7} = -aJ + J- + 2a J + - [a{N + 1) + 2b] J° - 2c J" + constant . (82) 

17 



For N > 0, there are N + 1 solutions which include excited states. Assuming <f>(z) = 
Y\^ =1 (z — Zi) in eq.(81), one obtains the Bethe ansatz equations which determine the roots 

ZiS 

z 2 

2az 2 - (2b + a)zi -2c-2a^ — = 0, i = l,...,N, (83) 

l^i Z i ~ Z l 

and the equation which gives the energy parameter in terms of the roots z^s 



N 1 



e 



2acJ2~ ■ (84) 



Each set of {z{\ determine a QES energy E with the corresponding polynomial 0. 

As an example, consider the j = 1/2 case with N — 1 and 0(z) =2 — 21. There are two 
solutions. From eq.(83), one sees that the root z 1 satisfies 

2az\ - (26 + a)z l - 2c = , (85) 

which gives two solutions 



± (2b + a)±^(2b + a) 2 + 16ac 
Zl = . (86) 



The corresponding energy parameters are 

1 



e ± = 2ac 

Zl 

a 



(2b + a)T J (2b + a) 2 + 16ac 



(87) 



For the parameters assumed here, the solution with root zj~ = — < gives the ground 
state, while that with root z{ > gives the first excited state. The superpotential W\ is 
constructed according to eq.(42) 

h(x) — z x 

= -ae~ ax + ce ax + + b . (88) 

1 + \z\ \e ax 

This gives the electric field 

= -aer ax + ce ax + 

1 + \z 1 \e 



pE™ (x) = -ae~ ax + ce ax + + b + k y . (89) 
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The ground state and the excited state have energy S 2 = m 2 + k 2 z and S 2 = m 2 + k 2 + 



In this paper we have shown that the Dirac-Pauli equation describing a neutral particle 
with anomalous magnetic moment in an external electric field is a QES system. Forms of 
electric field configurations permitting exact solutions, and QES states based on sl(2) algebra 
are classified in the spherical, cylindrical, and Cartesian coordinates. 

We have also demonstrated that the Dirac-Pauli equation with only electric fields are 
equivalent to generalized Dirac oscillators. Thus all the results obtained here are directly 
applicable to the later. Particularly, the exact solvability of the original Dirac oscillator is 
now recognized as one of the three exactly solvable cases of the Dirac-Pauli equation. 

It has been shown that the two-dimensional Dirac-Pauli equation is equivalent to a two- 
dimensional Dirac equation minimally coupled to a vector potential by a duality transforma- 
tion [14, 25]. However, as mentioned in Sect. I, the square of the Dirac Hamiltonian minimally 
coupled to a vector potential is proportional to the Pauli equation up to an additive constant. 
Hence the exact solvability discussed in [23] and the quasi-exactly solvability discussed in [12] 
of the Pauli equation can be straightforwardly applied to the two-dimensional Dirac-Pauli 
equation. 
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a J (26 + a) 2 + 16ac, respectively. 
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